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Abstract 

Special case of a Gibbsian facet process on a fixed window with a 
discrete orientation distribution and with increasing intensity of the un¬ 
derlying Poisson process is studied. All asymptotic joint moments for 
interaction U-statistics are calculated and using the method of moments 
the central limit theorem is derived. 


1 Introduction 

In the present paper we use methods developed in [9] to calculate all moments 
of Gibbsian {/-statistics of facets in a bounded window in arbitrary Euclidean 
dimension. These moments are used to derive the central limit theorem for 
such statistics. Central limit theorems for {/-statistics of Poisson processes were 
derived based on Malliavin calculus and the Stein method in [7]. 

Our calculations are based on the achievements in [T], where functionals of 
spatial point processes given by a density with respect to the Poisson process 
were investigated using the Fock space representation from [3]. This formula 
is applied to the product of a functional and the density and using a special 
class of functionals called {/-statistics closed formulas for mixed moments of 
functionals are obtained. In processes with densities the key characteristic is 
the correlation function j3] of arbitrary order which is dual to kernel function 
of the density as a function of the Poisson process. 

As in we call facets some compact subsets of hyperplanes with a given 
shape, size and orientation. Natural geometrical characteristics of the union of 
the facets, based on Hausdorff measure of the intersections of pairs, triplets, 
etc., of facets form {/-statistics. Building a parametric density from exponential 
family, the limitations for the space of parameters have to be given, so called 
submodels are investigated. In application of the moment formulas we are inter¬ 
ested in the limit behaviour when the intensity of the reference Poisson process 
tends to infinity. 
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We restrict ourselves to the facet model with finitely many orientations cor¬ 
responding to canonical vectors. In [9] basic asymptotic properties of the studied 
{/-statistics are derived. When the order of the submodel is not greater than 
the order of the observed {/-statistic then asymptoticaly the mean value of the 
{/-statistic vanishes. This leads to a degeneracy in the sense that some orienta¬ 
tions are missing. On the other hand when the order of the submodel is greater 
than the order of the observed U -statistic then the limit of correlation function 
is finite and nonzero and under selected standardization U -statistic tends almost 
surely to its non-zero expectation. By changing the standardization, however, 
we achieve a finite non-zero asymptotic variance. In the present paper we sim¬ 
plify the calculation of moments so that we are able to calculate any asymptotic 
moment. 


2 Central limit theorem 


Let Y = [0,6]'^ X {2b} x {ei,... ,ed} be a space of facets (facets are d — 1 
dimensional cubes) with three parts: set of facet centres, possible sizes of facet 
and possible orientations of hyperplane containing facet, i.e. we consider only 
facet with fixed size and orientation described only by elementary vectors. Facet 
in such case can be described as 

((zi,.. .,Zd),2b,ei) {{xi, ...,Xd)e M.‘^,xi = zf, la;^ - Zi\<b,i€ [d] \ {{}}, 

(zi,..., Zd) G [0, 6]'^, I G [d], [d] = {1,..., d} 

Moreover let (N, A/”) be a measurable space of integer-valued finite measures on 
Y, where Af is the smallest cr-algebra which makes the mappings x i—^ x(A) 
measurable for all Borel sets A CY and all x G N. We denote rja finite Poisson 
process of facets with intensity function aA on T, a > 1 in form 


aA(da;) = aA(d(z, r, (j))) 


= ax{z)dzS2b{r)-^Set{dcj)), 


where we have fixed the facet size, uniform distribution of the facet orientation 
and X : t K-i-. We also define interaction {/-statistics (using Hausdorff 

measure of order d — j) 

and the process fia with density 

/ d 

p(x) = Ca exp ^ r'jGj(x) 

\i=i 


with respect to ria, where a > 1, r'i is a real parameter and Cq = =- ^^ . 

Fullfilling of condition Vi < 0,* = 2,...,d assures that p(?7a) G fl 
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L?{Pn^)^ we will discuss necessity of these conditions later. We also use the 
notion of a submodel where Vj = 0, j ^ I and vi < 0. We will explore 
properties of such submodels with the order higher than 1, because in case 
we deal with Poisson process [S] . 


We can use a short expression for moment formulas using diagrams and 
partitions, see m, m- Let be the set of all partitions {<^ 2 } of [A^], where Ji 
are disjoint blocks and UJi = [k]. For k = ki + ■ ■ ■ + and blocks 

Ji = {j ■ ki -\ -h ki-i < j <ki-\ -h ki}, i = 1,..., TO, 

consider the partition tt = {Ji, 1 < i < to} and let H/d k,n Ilfc be the set 
of all partitions a € such that | J fl J' | < 1 for all J € tt and all J' G a. Here 
I J| is the cardinality of a block J G a. We will be referring to blocks of tt as to 
rows and we denote S{a) the number of pure singleton rows of partition a, i.e. 
number of J € tt, which for all J' G a, | J fl J'| = 1 ^ \ J'\ = 1. 


For a partition a G Hfei measurable functions fi : B ^ R,j = 

1,..., TO, we define the function —>• IR by replacing all variables 

of the tensor product that belong to the same block of cr by a new 

common variable, \a\ is the number of blocks in cr. We denote = 

where i repeats rui times for * = 1,..., s. It holds 0 =0,g>p. 
Now we state the main theorem of the paper. 

Theorem 1 Denote Gj{fx"a'^) = ^ l<j<d, 2<c<d, then 

2 

(Gi(Mi^)),...,Gd(AcM)) AZ, c = 2,...,d (1) 


as a tends to infinity, where Z ^ iV(0, E), E = {dij}ij=i, 
(c - 1) /c - 2\ /c - 2^ 


Ski = 


^k+i-i U_ W U _ 1 


Ikl, 


Ikl = 


2&, ei))W^-\di^fis^+k-i,2h, e.) n (si, 26, ei))x 
X x(si)dsi,..., x(sfc+i_i)dsfc+i_i, 


moreover 


G,(/rW) 0 , 

ai ^ di \ j ) 


c G {2, .. .,d},j > c, 
c G {2, .. .,d},j < c. 


where = /([o,6]d)j 2b, ei))x(si)dsi,..., x{sj)dsj. 
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Remark 1 Random variables GcifJ-^^),Gc+i{fJ-^^), ■ ■ ■ ,GdifJ-^^) are asymptot¬ 
ically degenerate, i.e. their expeetations tend to zero, thus eovariance of these 
variables Oki = 0,k > c,l G [d]. 

Remark 2 For random vector (Gi{r]a), ■ ■ ■, Gdipa)) we have similar results fS\/ 
with 0ki = dk+i~i 


Corollary 1 It holds 


■_ 1 



2,...,d,j <c 


as a tends to infinity, where Z ^ N{O,0jj). 


( 2 ) 


Lemma 1 It holds 

Pp{xi,..., Xp, 


Eexp{i^cGc{T]a U {xi, ■. ■ ,a;p})) 
Eexp(^cGc(j7a)) 


(d-ky. 


d\ ’ 


(3) 


as a tends to infinity, where Xi GY and k is number of facet orientations among 
{xi ,..., Xp} and c > 2. Moreover 


Pp{xi, ■ ■ ■, Xp, 


(d-fc)! 
(c—1 —fc)! 
d\ 

(c-l)! 


< Re"^“,R,S' > 0, 


where R and S do not depend on xi,... ,Xp 


Remark 3 The function Pp is called the correlation function of order p 

Proof: First consider submodel and facets Xi,... ,Xp with p < c distinct 
orientations and without loss of generality consider orientations ei,..., Cp, be¬ 
cause the distribution of orientations is uniform and the Pp{xi,... ,Xp, p^'^) does 
not change under rotation uniformly applied to all facets xi,... ,Xp. It holds 

Ppixi,. . .,Xp,p^a^) = 

_ Ylu=o ^ Iy” (^cGc{mi, ■.■,u„,Xi,..., Xp}) A"(d(ui,..., u„)) 

En=0 S /y" (I'cGclui, ..., Un}) A"(d(ui, . . . , Un)) 

We can obtain bounds for this expression by using the bounds for the volumes 
of intersection of facets b‘^~^ < < (26)'^“° as follows 


Z^n—0 n\ Z-/ni+...+nrf=n Vn 


) exp d, 


y^oo (^)" y^ 

Z_-/n=0 n! Z_-/ni + ...+nd 


< 


=n exp d, nG))) 

Pp{xi, . ■ . , Xp, p^^) < 


E" 0 '4f- E„.+,„«,.,. d. n<-))) 


y^°“ (-V) 

A^n—0 n! 


■E 


ni+...+nci. 


=n exp nG))) 


( 4 ) 
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where T = x('2)d2, Ui are the numbers of facets among ui,... ,Un with 

orientations e,, respectively, i = 1,... ,d and = (ni,..., rid). Furthermore 
denote 


Y[nj. 

FC[d] jeF 
c-P<\F\<c 
|Fu[g]|+p-g>c 


Specially d, is the total number of intersections of all c-tuples of the 

facets among ui,... ,Un and d, is the total number of intersections of 

all c-tuples of the facets among facets ui, cci,Xp. Then we substitute 

^ for a, extend the both fractions by and then we get in the case of 

the lower bound of (|1]) 


E CO 

ni =0 ■ * ■ Z— 


ai + .. 


nd—0 ni!...nd! 


■ exp d, - a(c - 1)) 


E oo ^-^co 

ni=0 ■ * ■ 2-^nd—' 


aj_-f • 


0 rai!...nd! 


■ exp (^Vcb‘^~‘^R‘^’^{0, d, — a{c — 1)) 


and in the case of the upper bound of (|3]) 


E OO ^-^oo 

ni—0 ■ * ■ A^Ud—O 


^i + . 


■ exp {vy-'^R^’P{p,d,n<^^'>) 


E OO V^OO 

m =0 ’ ’ * A^rid—' 


11 + .. 


Ud—O ni\...nd\ 


■ exp (z/c(26)‘^“‘^i?°’°(0, d, 


-q(c- 1)) 

-a(c- 1)) 


Using Lemma [2] we get the limit of the lower and upper bound in the same form 

(d —p)! 

■ For d > p > c we can get an upper bound in dH) just by using p = c, 

(C 7 I)! 

which tends to zero. 

Now consider more than one facet with the same orientation among xi,.. .Xp 
and with k < c distinct orientations, which are without loss of generality set to 
Cl,..., Cfe and maximum number of facets with the same orientation is P, then 
we can bound the correlation function 


/ aT \Ti 

E OO {- 3 -) 

n—0 n\ 


=n L nj exp {v,P’^{2bY-^R-’P{p, d, n(^))) 


ni-\-...-\-7id—Ti \ni,...,nd 


E“ 0 ^ {„z,j »P PM-R-ZO. d. nM>)) 

Pp (^1 ; • ■ • ; ; A^a ^ ) — 

Er.o ^ «P d. nW)) 


< 


A^n—0 n! A^ 


=n exp (fcP'^( 2 &)^-=PfO( 0, d, n(d-'>)) 


ni+...+nj. 

These bounds lead to expressions in the same form as in the case with dif¬ 
ferent orientations and thus we proceed in the same way and get the value of 

(d-k)! 

the limit For d > k > c we need only lower bound for the number of 

(P-l)! 

intersections in form d, which forms upper bound for the correla¬ 

tion function, which tends to zero. 
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Bounds for the numerator and denominator of the correlation function converge 
to their limits with at least exponential rate and we can also see that the upper 
bounds can be selected to depend only on the v, s and c , therefore they do not 
depend on currently selected facets xi,... ,Xp in the argument of correlation 
function. The rate of convergence can be extended to the whole fraction in [3] 


following way 









IA(a) 

-^1 


< i?2e 

-S2tl 

A{a) 

A 

A(a) 

^ l-fi2e-*'2“ 

- A 


— Sa 

B(a) 

B 

^ B 


^ (1 - i?2e-‘52“).B 

— 


A{a) 

A ^ 

A(a) 

l-|-fl2e-“2“ 

■A 

^ig-Sia _ AR2e~^^°- 


-Sa 

B{a) 

B 

B 


(1 -b i?2e-^2“)B 

J LC' 
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where A{a) is the value of numerator and B{a) is the value of denominator on 
the left side in ([3|) ,respectively A, B is the limit of the numerator, denominator 
on the right side in ([3]). 


Lemma 2 


E-E 


a 


ni + ...+nd 


ni! ...nA 


exp {vR‘^’P{p, d, - a(c - 1)) 


m—O nfi—O 

where v<0, p<c<d, c>2 and i?, S' > 0. 


{d-p)\ 
(c- 1 - p)! 




Proof: We denote 

°° °° „ni + ...+ns , . 

I{a,c,t,s) = ^ —j-^ exp s, ) - a(c - l)j . 

rti=0 ns=0 

Firstly we calculate the values of the limit by calculating the sum over 
(ni > 0 A ... A > 0) to show that this value tends to zero as a tends to infin¬ 
ity. We show this only for p = 0 because for p > 0, we get upper bound using 
p = 0 and the sum is non-negative. In the following we use Chernoff bound for 
tail probabilities of Poisson distribution 


m 


E 

1=0 



(es)™ 

to"* 


,m < s. 


1. First we consider that all the summing variables are between 0 and 
((a^/^ > m > 0) A ... A (a^/^ > Ud > 0)) 


ni>0 nd>0 


, 2/3 


, 2/3 


E-E 


a 


ni + ...+nd 


m! ...udl 


exp (—a(c — 1)) < 


(ea)“ 


ni>0 nd>0 

where we used d-times the Chernoff bound. 


(a2/3) 




g-a(c-l) ^ Q 
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2. Now consider that one of the summing variables is greater than 
without loss of generality we select nj, 

> ni > 0) A ... A > n^-i > 0) A {ud > 

aV3 „2/3 


E-- E E 


°° ^ni + ...+nd 


m'. ...nd\ 


ni>0 >0 n(i>a3/3 

2/3 2/3 


exp d, — a(c — 1)^ < 


E- E E 

ni>0 nd_i>0nd=0 

E- E 


a 


ni + ...+nd 


ni!.. .rid! 

ni + ...+nd_i 


ni>0 nd_i>0 


ml.. .nd-i! 


exp {vud — a{c — 1)) = 


exp (ae'' — a{c — 1)) < 




2/3 


exp {ae'' — a{c — 1)) ^ 0 


^(a2/3)(d-l)a3/3 

because e'' — {c — 1) < 0. 

3. When at least two of the summing variables are greater than without 
loss of generality select Ud-i and nd, then we have 

£ j 2/3 ^ 2/3 ^ ^ +nd 


£2. ^ni + .-.+Kd 


ni>0 ^d- 2 >0 1 >a3/3 nrf>a3/3 

ill--- E E E 

ni>0 nd-2>0nrf_i>a2/3„j>a2/3 

< exp + a{d + 1 — c)^ —>■ 0. 


exp — a(c — 1)^ 


4. The same applies to the case, where more than two variables are greater 
than because we are able to find terms with higher power of a in the 
exponential. 

Therefore we need only to examine the remaining terms, where at least one of 
the variables is equal to zero, thus we replace I{a,c,p,d) by d sums, where one 
variable is set to zero 


/(a, c,p, d) « p/(a, c,p — 1, d — 1) + (d — p)I{a, c,p, d — 1), (5) 

where w is the equality after omiting the summands, which tend to zero on the 
left side, I{a,c,p— 1, d — 1) is the sum after setting to zero one of the variables 
rii, ..., Up, I{a, c,p,d — 1) is the sum after setting to zero one of the variables 
rip+i ,... ,nd and the multiplying numbers are the counts of possible selections 
of these variables. It can be shown that 

lim I{a,c,t,c — 1) = 0, t < p, 

a—^oo 

= 1, t = p. 
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Because the series on the right side of ([S]) are in the same form as the original 
one and we can again sum only over the indices, where at least one is equal to 
zero, thus we repeat the {d — c + l)-times step in @ and we get 

d—c-\-l 

I{a,c,p,d)K. ^ Cj/(a,c,p-j,c-1), (6) 

j=o 

where Cj € N. All summands tend to zero with one exception of 
co/(a, c^p,c— 1) with 


Co 


(d-p)! 

(c-l-p)! 


,c>p, 


0,c = p, 


which is the number of possible selections of variables set to zero from rip+i ,... ,nd 
in d — c + 1 steps. The overall speed of convergence is implied by the conver¬ 
gence speed of every part of the sum, which converges to its limit at least at 
exponential rate. □ 


Lemma 3 It holds 




{^j=0 ^ (mi, . . . , M|,t|)P|,t| (mi, ■ • ■ , (d(wi, . . . , M|^|))- 




(Ml,..., W|.|)Al'^l(d(ui,..., W|.|)) 


< Re 


— Sa 


where a G S,R> 0, Yc-i = [0, x {26} x {ei,..., Cc-ij is space 

of facets with d — c-I- 1 orientations missing (we can select any orientations), 
mj G N° and ..., Uj) = . 


Remark 4 In space Yc-i we still use intensity A, even when the orientations 
are restricted to {ei,..., Cc-ij and the measure is non-zero on orientations 
(ci,..., Cd). 


Proof: 

The limit of correlation function depends only on the number k of the distinct 
orientations among the facets (iti,... ,U|ct|), then correlation function tends to 





(d-fc)! 

and thus we can write 

(<=-!)! 


/yki 


^ (ui,...,(ui,...,ni^i,(d(ni,...,«|<,|)) 


E 


E 

/c=l 


kj J(Y\-\\u] 
xp|^|(ui,... ,U|„|,/ri'=^)Al'"l(d('Ui,... ,W|<,|)) < 

{d-k)\ 


d 

E 


(®,to (Ui,... ,n|.|)^^f^Al-l(d(ni,..., n|.|))+ 


(c-l)! 


fe=l 


/(Fi-i; 


[fe] 


{d-k)\ 


Al‘^l(d(Mi,...,M|^|)) < 


E 

k^l 

d 

E 

fc=l 


c — 1 




(®to (^1’ • ■ • .«kl)A'"'(dK,..., W|.|))+ 

i)w ^ ^ 

(®,to (^^1,..., U|^|)Al"l (d(z.i,..., U|^|))i?e-^“ 


where is subspace of where facets ui,...,u\cr\ use orientations 

ei,..., Cfe (each orientation is used at least by one of the fasets), (^) is the number 
of possible selections of orientations used. We have an upper bound for the 
expression in the absolute value and we can get a lower bound in the same way. 
□ 


Remark 5 This means that we do not need the correlation function in the 
calculations and we can omit it from the formula, if we also omit d — c + 1 
possible orientations. 

Proof of Theorem 
It holds [I] 


J- JYo 


c-l 


c-l 




El 


i=i 




Inij 


(«i, ... ,«|<,|)p|.|(ui,... E^)Al'^l(d(«i, 


(7) 

( 8 ) 

■ I 'a|(T|))- 
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We can also get a relation for joint moments of centered random variables 


c-l 


E 




i=i 


j=i 


ii—0 ic-i—0 


mi 

h 


rric-i 

1 




xE 





nE(G,(Mi^)))*\ 


(9) 


where = M. 

Firstly we calculate expectations of the U-statistics using Lemma [3] and 





A^(d(ui,...,uj)) = 


where j^) is number of possibilities how to select unique j orientations from 
c — 1 and jl is number of possible allocation of them on the j positions, P 
is number of all possibilities how to select orientations from d possibilities. Ij 
is value of integrated Hausdorff measure over the space of centres, which does 
not depend on the currently selected orientations, they only need to be distinct, 
otherwise the measure would be 0, i.e. only non-paralel facets intersect and thus 
have non-zero volume of the intersection. We can already see that expectation 
of {/-statistics of order higher or equal than c are zero, therefore they converge 
to zero in and thus we only need to investigate the {/-statistics of the order 
lower than c. 

Secondly we calculate all joint moments. To do this we need to first use formula 
([3]) and Lemma [3] - we use the limit values of correlation function, which we 
justify later. To describe the relation between original formula and the formula 
with correlation function replaced by its limit value we use ~ 


' c-l 


'c-l 


c-l 




c-l 


U=1 




:n H"-^(nLiZiOA^(d(ui,...,u,)) 

j=i / 


E 






(m-i - ii - ic-l) 




(^^1, ■ ■ ■, «kl)A'"' (d(ui,..., u|<,|)) 


We are interested only in terms with power higher than or equal to M, 
because the other terms will tend to zero with increasing a, i.e. partitions 
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fullfilling condition \a\> M — YTj=i do not have to examine odd 

moments, i.e. those with odd, because there is not any summand with 

the power of a matching M in the divisor, thus they can be only zero or infinite, 
therefore if we prove that all even moments tend to some finite value, then all 
odd moments are equal to zero. 

Select s = (si,..., Sc-i), so that > Si > 0, i G [c — 1], 3j G [c — l],rnj > sj, 
choose any partition Us G fullfilling conditions |(Ts| > 

S{(7s) = 0, i.e. each block of tt is connected to any other block of tt by some block 
of s. Then for t = (ti,..., tc-i),'rni > U > Si,i G [c — 1], 3j G [c — 1], > Sj 

there are partitions at G 11* which have only additional singleton rows 

compared to CTs, S{at) = J2i=i ti ~ Si = |crt| — |(Ts| and it holds 





[ H‘'-dnyiuE^(d(wi,.. 

.,Uj)) 


J 


because we can separate the singleton rows corresponding to the functions 
in tensor product, which can be integrated separately, because they do not have 
any common variables with the other functions in the tensor product and the 
integral is equal to the expectation of {/-statistic. We can see that all summands 
corresponding to any of the partitions at in the evaluation of (jH]) contain common 
term 


0 = (®5=i {uu . ■., ..., 

n ( (nLi«.) A^ (d(ui,..., u,)) j 


i=i 

and then we sum over all such partitions at 

mi mc — i 

0 E-- E 


2l_Si lc-l—Sc-1 


nil 


rric-i \ / ti 

ic-i J Vsi 


Sc-lJ 


0(-l) 

mi —Si Sc-i 


—1' 


TOi 

Si 


rric-i 

Sc-l 


mi — & ! ''‘'c — i'=c—1/ \ / \ 


where we use Binomial theorem for summing with necessary condition ^ 

YTjZ^ and (y) are original coefficients from formula (0 and (*^.) is the num¬ 
ber of options how to select additional singleton rows. 
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Therefore all partitions with any singleton rows or containted within Ilf ,,-1) ® < 
m = (mi,... ,mc-i) cancel each other out. But we calculated with the limit 
values of correlation functions and the integrals are multiplied by a in polyno¬ 
mial, thus we have to deal with speed of convergence. We have the limit in a 
form 


N „ ^ f 

X! / H{x.)Qi{x,a)A{x.) < a'' X! / ^^(x)|£'i(x,a) - ft(x)|A(x) 

i=0 1=0 


< 


, ^ r 

< a^Re-^^ / ^(x)A(x) 


i=0 


lx 


where H >0 , Qi > 0,|pi(x, a) — Pi(x)| < Re , R, S > 0, 

Ix -^(x)ft(x)A(x) = 0, the sum is over all partitions, which are the same 
after leaving out all singleton rows and the measure represents product measure 
of A on space X = YP for some p and K is created by multiplying Hausdorff 
measures and Qi by multiplying correlation functions. 

Now we are left out only with partitions cr, which do not contain any pure single- 
ton rows, are contained in These partitions have each row connected 

exactly to one another row by one block of two elements in a (|ct| = M) and 
therefore, if we omit all the mentioned partitions, then 



c-l 


-M 


T. T. T. n 

=1 <T6nK,Je<T:| J|=2 J'={bl,&2}6o' 


0-1 


i(d(a;i,...,x^(f,j)+^(b,)_i)),T(s) = maxj-g[c_i] < ^m* < s = ^ 


rrii 






where we sum first over all possible selections of common elements among the 
partitions and then over all possible pairings of partition rows, we also divide 
integral into several parts, where each part consists only of elements which are 
in the same block of a partition. Function r connects row of partition to its 
length. It holds 


(c- l)(r(i„) - 1 )!(t( 6.) - l)!/„6,„(»,|(,yA.)(,yf_.) 


d'r(6i)-|-T(62)-l 


( 10 ) 
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where c — 1 is number of possibilities how to select the one common facet ori¬ 
entation, (-r(6~)^i) is number of possibilities how to select distinct re¬ 

maining orientations of the rest of the facets in the first and the second function 
in integrand and (t(&i) — 1)!,(t(&2) — 1)! are numbers of their possible positions 
among facets in the Hausdorff measure, jg number of all possible 

orientations from all d possibilities (even non-distinct ones) and ^T(bi)r(&2) 
integral over facets with fixed orientations over the space of the facet centres. 
Then using (fTUl) 

E {u-zl n;=} E ^ 


(5) 


M 


E 


E 


c-l 

E 


n 


<TGn_K-,jg(T:|j|=2 

n (-(0.) - 1) 

c-2 \/ c-2 \ 

^(^1) - V V^(^2) - V ■ 


(5)" i: n 


<7eUK,Je<y:\J\=2 J={bi,b2}eo- 


If we express the covariance of any two variables in the same form by selecting 
and and then by calculating 

a*-st-i(c _ l)/„. (c -2\{c- 2\ 

dFJ^^ 0 - 0 0 - V 


we can see that the distribution of statistics has the property of normal distri¬ 
bution, i.e. joint moments of centered variables are equal to sum over all pairs 
of unordered random variables (random variables with higher power are used as 
several distinct multiplied random variables) and this implies the central limit 
theorem, because normal distribution is defined by its moments [2]. 

There is only one remaining statement to prove 


ai ^ di \ j ) 


cG {2,...,d},j < c. 


the first moment of the random variable on the left side is equal to right side 
and the variance tends to zero as can be seen from the central limit theorem. □ 


Remark 6 Consider process pa with density in more general form p(x) = 
Ca exp ( b'iGi(x) j. Assume there is c > 2, i^c > 0 o.nd select such mini- 
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mal c. 


Eexp '^VjGj{r]a) = 


0=1 


^n^—aT r I ^ \ 

n=0 ”■ \j=l I 




°° ^ni+...+nd 


ni=0 nd=0 


^^-“"E-E 

ni=0 ric—O 


ml ...udl 


gni+...+ra<; 

ml ■■ ■ ricl 


E 11’ 

j=l {ii,...,ij}c[d] 1=1 

“pIE-; E ri' 

0-1 {ii....,ij}C[c] i=l 


S»-“"E7^p’<P Ep; 


^ (n!) 

n=0 1 ^ 


0=1 


( 11 ) 


where v'- Vj > 0 and we firstly set the last d — c summing variables to zero and 
then summed only over the summands, where all of the summing variables have 
the same value. It can be proven (e.g. by using ratio test), that the sum in 
m is divergent and therefore p(r]a) (f L^{Pr)fi) in this case. On the other hand 
non-positivity of parameters v implies p(j]a) G L^{Pv ) riL^(P,, ) as shown in 
which finally leads to vi <Q,l >2 viVa) G L^iPva) ^ P'^iPva)- 


Remark 7 Consider process p-a with density p(x) = Ca exp t'iGi(x) j, 

E 0, / > 2. Assume there is c > 2,^ < ^ and select minimal such c. Then 
using similar techniques as in proof of LemmoQ] and Lemma\^ we can show that 

lim Pp{xi,.. .,Xp,pa) = lim Pp{xi,.. .,Xp,p‘fi'>), 


which leads to the same asymptotic distribution of statistics (Gi{pa), ■ ■ ■, Gd{pa)) 
as {Gi{pi^^),...,Gd{pi^^)). 
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